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ABSTRACT 

The primary focus of this paper is on the particle acceleration mechanism in solar coronal 
3D reconnection null-point regions. Starting from a potential field extrapolation of a SOHO 
magnetogram taken on 2002 November 16, we first performed MHD simulations with horizontal 
motions observed by SOHO applied to the photospheric boundary of the computational box. 
After a build-up of electric current in the fan-plane of the null-point, a sub-section of the evolved 
MHD data was used as initial and boundary conditions for a kinetic particle-in-cell model of 
the plasma. We find that sub-relativistic electron acceleration is mainly driven by a systematic 
electric field in the current sheet. A non-thermal population of electrons with a power-law 
distribution in energy forms, featuring a power-law index of about -1.75. This work provides a 
first step towards bridging the gap between macroscopic scales on the order of hundreds of Mm 
and kinetic scales on the order of cm in the solar corona, and explains how to achieve such a 
cross-scale coupling by utilizing either physical modifications or (equivalent) modifications of the 
constants of nature. With their exceptionally high resolution — up to 135 billion particles and 
3.5 billion grid cells of size 17.5 km — these simulations offer a new opportunity to study particle 
acceleration in solar-like settings. 

Subject headings: Corona — DC electric field — particle acceleration — Sun 



1. Introduction 

During solar flares an enormous amount of 
energy is released, in particular in the form of 
highly energetic non-thermal electrons. It is gen- 
erally accepted that the underlying release mech- 
anism is magnetic reconnection. In connection 
with a reconnecting current she et, strong large- 
scale electric fields can build up (jLiu et al. f (l2009h 
and references therein), leading to direct accel- 
eration of particles beyond thermal velocities, 
while fluctuating electric fields created by re- 
connection and other dynami c events can lead 
to stochastic acceleration ( e.g. Miller et"aLlll996l 
19971: IPetrosian et al.|[2006l ). 

During the last decade high resolution ob- 
servations of solar magnetic fields from several 



solar space missions and ground-based obser- 
vations, such as YOHKOH, SOHO, TRACE, 



LouECODe et al.l2005l Ko et al.ll2003; Milhgan et al. 


20061: Priest & Schriiver' 


1999': 'Sui et al. 2004; 


Aulanier et al. 2007: Kumar et al. 2011i: 


Jess et al. 


2008) have brought new insights and substantial 



support for existing solar flare and CME models. 

Additionally, new three-dimensional coronal 
mag netic reconnection and acceleration mo dels 
(e.g. IPriest fc Titovlll996t IVlahos et all 12004 ) as 



well as the rapid expansion of c omputing resources 



for large scale siinulatio ns ( Isobe et al] 2007t 



Galsgaa rd fc PontinI 1201 It iToriumi &: Yokovama 



12012,) have opened up a new chapter in the un- 
derstanding of the formation of current sheets and 
particle acceleration sites in three dimensional re- 
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connection regions. But what is lacking is an 
interconnected understanding of how microscopic 
plasma physics scales interact and exchange in- 
formation with macroscopic large-scale magne- 
tohydrodynamics (MHD) scales in the solar at- 
mosphere. MHD scales provide the environment 
for particle acceleration to happen, but it is un- 
clear how the complex non-linear feedback from 
much smaller scales onto the overall behavior of 
the plasma above the solar surface is handled by 
nature. 

Most studies have made use of the fact that 
the temporal evolution of the large-scale mag- 
netic field in the solar atmosphere can to a 
first approximation be described by compress- 
ible MHD. But fluid approaches are limited to 
thermal particle distributions, which is not suffi- 
cient to describe the kinetic aspects of magnetic 
reconnection that convert magnetic field energy 
into particle kinetic energy. A proper descrip- 
tion of such processes requires taking into ac- 
count the back reaction of kinetic processes on 
the large scale dynamics. N evertheless, test parti- 
cle MHD simulations, ( e.g. Turkmani et al. 20051 
2006|: iBrowning et all I2OIOI: iDalla fc Browning 



2005! I2OO8I: iRosdahl fc GalsgaardI l2010l ). gfrc a 
good idea of the overall acceleration region frame- 
work. In such simulations MHD fields evolve in- 
dependently of the motion of the test particles. 
There is no immediate backreaction from the ac- 
celerated particles onto the fields, which is po- 
tentially a serious limitation, since especially the 
changes of the electric field that would be induced 
by the accelerated particles can be large com- 
pared to the backgr ound field induced b y mag- 
netic reconnection ( Siverskv fc Zharkoval [2009). 
The lack of feedback can lead to an exaggerated 
particle acceleration, as has been noted by for ex- 
ample [Rosdahl & Galsgaard (2010), since there 
is no limitation for the energy gain of particles. 
Furthermore, kinetic instabilities can be of impor- 
tance for the fast reconnection onset in solar flares, 
and more generally for the evolution of the cur- 



rent sheets in reconnection regions (jBarta et al 



There is therefore a need for realistic self- 
consistent kinetic simulations to examine micro 
scale processes in plasmas and to be able to prop- 
erly take into account backreactions from the par- 
ticles to the flelds. 

The main challenge for kinetic simulations on 



scales of solar events is the enormous dynamic 
range involved. Explicit particle-in-cell (PIC) 
simulations have to resolve characteristic kinetic 
scales and are restricted to very small physical 
sizes, due to limitations that arise from condi- 
tions for code stability as well as from resolution 
criteria. These computational restrictions have 
so far prevented investigations of the coupling of 
kinetic to MHD scales using kinetic simulations, 
and have thus prevented self-consistent modeling 
of particle acceleration in solar flares. There has 
m ainly been one interlocked m odel attempt e.g. 
bv ISugivama fc Kusano ( 20071 ). but the complex- 
ity of the problem limits the applicability of such 
attempts significantly. In the present study we in- 
troduce a new method to meet the challenge of 
multi-hierarchy simulations and use the method 
to investigate particle acceleration mechanisms in 
a solar reconnection event using ultra large scale 
kinetic modeling. 

In Section[2]wc describe the numerical methods 
and their implementation, while in Section[3]we in- 
troduce the experimental setup and list the most 
important simulations we performed, indicating 
their relative roles and importance. In Section 3] 
we present and discuss the results, and in Section[S] 
we summarize the results, present our conclusions, 
and give an outlook onto future work. 

2. Methods 

We perform PIC simulations using the Photon - 
Plasma code (iHaugboUd 120051: iHededall l2005l) . 
which solves the Maxwell equations 



dB 

'dt 
9E 



= V X E 



eoAio-^ = V X B - ^oJ, 



(1) 
(2) 



together with the relativistic equation of motion 
for charged particles 

^cf(7v) ^ ^(E + vxB) (3) 
at 



on a staggered Yee lattice ( Yed 1966 ). We use SI 
units, scaled so the unit of length is 1 km, the unit 
of time is 0.1 s, and the unit of density is 10 kg 
m~^. 

The Lorentz force is computed by interpola- 
tion of the electromagnetic fields E and B from 
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the mesh to the particle positions, employing a 
cubic scheme using the 64 nearest mesh points. 
The code integrates the trajectories of protons and 
electrons moving in t he electrom agnetic field with 
a Vay particle mover (|Vavll2008l ). The charge den- 
sity p is determined by weighted averaging of the 
particles to the mesh points using the same cubic 
scheme as for the field to particle interpolation, to 
avoid self forces. The currents are found using a 
new 6'^ order versi on of the Esirkep ov charge con- 



servation method (|Esirkepovll2001l ) that is consis- 
tent with the field solveiQ- The field equations are 
solved on the mesh using an implicit 2"*^ order in 
time and 6*'' order in space method. Because the 
solver is charge conserving and the fields properly 
staggered, Gauss' law is obeyed and the B field is 
kept divergence free to numerical precision. 

The boundaries of the domain are fixed for the 
magnetic fields and open for particles. Particles 
can escape and new particles are added to the box 
from 'ghost cells' outside the physical boundaries, 
where the conditions are specified from values in 
the MHD snapshot. In the current short duration 
kinetic simulations the field values in the bound- 
aries are held fixed. In longer duration simulations 
they could be made time dependent by performing 
interpol ations in time b e tween MHD snapshots. 

(l2012l) . 



See also Haugbolle et al 



Below we introduce for the first time a modifi- 
cation of the elementary charge q, in addition to 
and analogous to the well known speed of light 
(c) modification th at has previously been used 
in many cases (e.g. IPrake et al. 2006 ). Since all 
the micro-scales (gyro radii, skin depths, Debye 
lengths) are inversely proportional to the charge 
per particle, one can increase the micro-scales un- 
til they are resolvable on macroscopic scales by 
decreasing the charge per particle sufficiently. 

While changing the ratio of micro- to macro- 
scales by a large amount may appear to be a very 
drastic approach, the method can be defended on 
both qualitative and quantitative grounds: From 
a qualitative point of view it may be argued that 
as long as one retains a proper ordering of non- 
dimensional parameters, as discussed in more de- 
tail below, one should expect to see essentially the 



same qualitative behavior, albeit with (possibly 
large) differences in quantitative aspects. From 
a quantitative point of view one can indeed at- 
tempt to predict how these quantitative aspects 
depend on the modifications of q and c, to be able 
to extrapolate — at least to order of magnitude — 
to values that would be typical in the unmodified 
system. 

Adopting this approach enables us to perform 
explicit PIC simulations of large scale plasmas. 
With the parameter values used here we resolve 
the electron skin depth 5e with at least 3.8 grid 
cells and the Debye length with at least 0.3 grid 
cells. For the time stepping we use a Courant con- 
dition of 0.4, considering the light crossing time in 
a grid cell as well as the local plasma frequency. 

3. Simulations 

In order to investigate the particle accelera- 
tion mechanism around a 3D reconnection region, 
we started out with a Fourier transform potential 
extrapolation of a Solar and Heliospheric Obser- 
vatory/Michelson Doppler Imager magnetogram 
from 2002 November 16 at 06:27:00 UT, 8 hour 
prior to a C-fiare occurrence in the AR10191 
active region . Thi s is similar to the setup by 
Masson et all (|2009l) . As there is no vector magne- 
togram available for this event, a non-linear force- 
free extrapolation has not been feasible. The po- 
tential field magnetic configuration arising from 
the extrapolation showed mainly two connectivity 
regions, separated by a dome-shaped 'fan-surface' 
(jCraig et al.l 119971 ) . each including a spine struc- 
ture, which is the symmetry line intersecting the 
fan at the magnetic null-point. 

We define the inner spine as the spine inter- 
secting with the solar surface inside the dome-like 
surface delimited by the fan-plane, while the outer 
spine reaches up into the corona before return- 
ing to the photosphere several tens of Mm away 
from the null-point (cf. Figure [1]). The magnetic 
field lines in the fan-plane meet at the magnetic 
null-point, which is located at about 4 Mm above 
the photospheric boundaryQ The dome-like fan- 
plane is a typical feature of a parasitic polarity 



^The original Esirkepov method, just like most other charge 
conserving schemes in the literature, is only consistent with 
a 2""^ order field solver. 



^For computational reasons we tapered off the vertical mag- 
netic field towards the boundaries of our horizontally peri- 
odic MHD model, and the hei ght of the n ull-point is there- 
fore different from the case in iMasson et al.l (20091) . 
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magnetic null-point topology originating when a 
vertical dipole field emerges into a magnetic field 
configuration of opposite magnetic polarity. For 
an elaborated descriptio n of the observation s of 
this solar flare event, see Masson et al. I (I2OO9I) . 

We previously conducted fully 3D resistive 
compressible M HD simulations using the S taq- 
ger MHD code (|Nordlund fc Galsgaardlll997t ) on 
a box size of 175 x 100 x 62 Mm and a stretched 
896 X 512 X 320 grid with a minimum cell size 
of 80 km in a reg ion around the null-point 
( Baumann et al.l 12012 ). In most cases the initial 
particle density and temperature were set to the 
constant values of 6.8 x 10 cm" -3 and 5 x 10^ K, 
neglec ting gravity, analogously to .Masson et al. 
()2009l) . The boundaries were chosen to be closed 



at the bottom and top, and peri o dic on the sides 



Assuming, as did Masson et al. ( 2009[ ). that the 



strong horizontal motions observed by SOHO dur- 
ing the considered time interval was the driver of 
the reconnection event, we simulated the plasma 
motion caused by these motions by using essen- 
tially the same elliptically shaped driver at the 
ph otospheric boundary of the box as was used 
bv iMasson et al. ( 2009t ). applying velocity ampli- 
tudes ranging from 3.33 to 20km s~^ and a hor- 
izontal shape function schematically representing 
the observed photospheric motion. The enhanced 
driving speed, compared to observations, led to a 
significant reduction of the simulation time, while 
it was still sub-Alfvenic and thus left the overlying 
plasma and magnetic field enough time to adapt 
to the displacement. 

The line-tied photospheric boundary motions 
indirectly reshuffled the fan-spine geometry at its 
foot-points by a pressure increase onto the fan- 
plane, which caused a relative displacement of 
the magnetic field lines inside and outside the fan 
plane of the magnetic null point, respectively, re- 
sulting in the formation of a current sheet in the 
fan-plane. The displacement of field lines also 
causes a relative displacement of the two spines, 
which then leads to magnetic reconnection at the 
null-point and a growing resistive electric field in 
the fan-plane. A detailed description of the MHD 
simulations and t heir results are publi shed in a 
companion paper (jBaumann et al.ll2012l ). 

For the present study we used an MHD sim- 
ulation with an initial particle density of 6.8 x 
10^^ cm~^, a temperature of 5 x 10^ K and an ap- 



plied driving speed of 20 km s~^ . After 240 seconds 
of simulated solar time the null-point area showed 
clearly enhanced current densities parallel to the 
magnetic field, indicating that a dissipative pro- 
cess was taking place. We choose this as the start- 
ing point for the 3D relativistic PIC simulation. 
To minimize computational constraints due to the 
plasma frequency we rescale the density from the 
essentia lly chromospheric va lue of 6.8 x 10^^ cm 



used bv lMasson et al. ( 2009[) partway towards val- 
ues more characteristic of an active region corona, 
so that it becomes instead 1.28 x 10^^ cm""^. 

Figures [T]i and[T}D illustrate the initial PIC sim- 
ulation setup, by showing the chosen cut-out of a 
snapshot from the MHD simulation, as seen along 
the X and y axis respectively. The outer spine 
extends to the right of Figure [Dd, while the fan- 
surface spreads out over the rest of the area. Due 
to t he initial difference in the fan-plane eigenval- 
ues (iM asson et al.|[2009[l a slight asymmetry in the 
fan-surface is noticeable. The inner spine can be 
recognized inside the volume spanned by the fan- 
plane field lines in Figure [TJi. 

For the kinetic simulations we use the Photon- 
Plasma code on a uniform grid. A cut-out 
of the MHD simulation box with dimensions 
44x25x16 Mm is chosen for a set of simula- 
tions, covering cell sizes from 70 km down to 
17.5 km corresponding to uniform grids with up 
to 2518x1438x923 ceUs. The simulations are gen- 
erally performed with 20 particles per species per 
cell, covering up to 30 solar seconds. We neglect 
gravitation, since the aim of this study is mainly 
to assess the electron acceleration mechanism, for 
which the influence of the gravitational force is 
negligible. The initial ion velocities consist of two 
components; a random thermal velocity drawn 
from a Maxwellian distribution plus the bulk ve- 
locity from the MHD simulation. For the electron 
velocities, we use the sum of a random thermal 
velocity drawn from a Maxwellian distribution, 
the bulk velocity and the velocity due to the ini- 
tial electric current. The initial electric field is 
simply taken as the convective electric field from 
the MHD simulation, which will be self-corrected 
by the system after the first few time-steps of the 
kinetic simulation, by adapting the electric field as 
well as the electric current to balance the system. 
An overview of the relevant simulations may be 
found in Table [H 
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3.1. Modifications to the System 

It is common in MHD and PIC simulations to 
convert physical units into numerically more con- 
venient code units. Such a rescaling leaves the 
simulated physical situation completely unaffected 
and all parameters can easily be rescaled back to 
physical units at any time. Because they leave the 
ratio of scales unaffected these kinds of conversions 
between units of measurement are not sufficient to 
allow explicit PIC codes to address the multi-scale 
issues discussed below. 

Magnetic reconnection is a typical example of 
multi-scale physics, where microscopic scales inter- 
act with and couple to macroscopic scales. While 
kinetic scales in the solar corona are on the order 




Fig. 1. — Initial state of the PIC code simula- 
tion; a cut-out of size 44 x 25 x 1 6 Mm of an evolved 
MHD simulation snapshot, a) Seen along the x- 
axis and b) along the y-axis. Shown are magnetic 
field lines (purple), the maximum current density 
at the lower boundary (blue) and the electric cur- 
rent density as a semitransparent direct volume 
rendering (yellow for lowest value, orange for high- 
est value). 



of mm, coronal structures such as the ones investi- 
gated here have scales on the order of tens of Mm. 
This range of spatial scales is impossible to cover 
by explicit kinetic simulations — now as well as in 
the foreseeable future. Explicit PIC codes are sub- 
ject to a number of numerical stability constraints. 
For the code that we use here, one of them is that 
the Debye length, defined as Ad = a/ e^ksT /Uf-q^, 
where fcs is the Boltzmann constant, T the tem- 
perature, eo the vacuum permittivity, the elec- 
tron density and qe the electron charge, should be 
at least on the order of a fraction of a grid cell As: 

0.3 As < Ad. (4) 

In addition there are also stability constraints on 
the time step set by the speed of light and the 
plasma frequency 

At < As/c, (5) 

At < 2/wp. (6) 

In summary, code stability requires that the time 
step needs to approximately resolve the light-wave 
and plasma wave propagation and that the grid 
spacing should not severely under-resolve the elec- 
tron Debye length. When these conditions are 
marginally fulfilled in the global simulation do- 
main they are typically (except for the speed of 
light condition) fulfilled with a good margin in 
most computational cells, which helps to improve 
the accuracy of the numerical solution. 

Additionally we need to fulfill the following ap- 
proximate equality: 

Sskin ~ a few • As. (7) 

This requirement arises to ensure a reasonably 
faithful representation of the plasma, as employ- 
ing too few grid cells per skin depth suppresses 
sub-skin depth plasma behavior, while, for a fixed 
number of grid cells, having too many grid cells 
per skin depth results in a model where the kinetic 
scales and the MHD scale are not sufficiently sepa- 
rated, and consequently the MHD behavior is lost. 
Given these constraints, we would require approxi- 
mately one septillion (1 x 10^'*) grid cells for an un- 
modified simulation to simultaneously resolve the 
microscopic and macroscopic scales in the prob- 
lem. This fundamental constraint notwithstand- 
ing, bridging the large scale difference is in fact 
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approachable. One can either change some of 
the physical properties (e.g. the magnetic field 
strength or the plasma density) or, equivalently, 
modify some of the constants of nature. As long 
as such modifications of the physical properties of 
the system do not change the large scale behavior, 
they may be helpful by decreasing the gap between 
microscopic and macroscopic scales. 

We first introduce two possible sets of mod- 
ifications of the physical properties which con- 
duce to decrease the large gap between the micro 
scales and the macro scales, using three modifica- 
tion parameters, which are closely related to each 
other, so that only two of them are so to speak 
free parameters, while the third one is a result 
thereof. One of the reasons for this is an MHD 
variable constraint; we wish to maintain the rel- 
ative importance of magnetic and fluid pressure, 
i.e., we wish to keep the plasma beta parameter 
unchanged in each type of modification. The goal 
is to change micro scales, which cause the pre- 
viously mentioned numerical constraints, i.e. in- 
creasing the Debye length A d , decreasing the span 
of velocities (equivalent to decreasing the speed of 
light), and decreasing the plasma frequency. In 
modification of type A we change the magnetic 
field strength by a factor 6 < 1 , implying a change 
in the density to keep the pressure balance. This 
increases the length scales (A^, the gyro radius Vg 
and the skin depth 6skin) and at the same time re- 
duces the frequencies {ujp and the gyro frequency 
ujg) by the same factor, while leaving the temper- 
ature and hence the sound speed Vs and Alfven 
speed VA of the system unchanged. 

Modifications of type B address the macro- 
scopic speed ratios. We start by changing the 
temperature by a factor which then leads to 
an increase in the acceleration of gravity g by the 
same factor. The consequences of this modifica- 
tion are that speeds increase by a factor t, while 
the length and frequency ratios are forced to be 
changed to achieve these required modifications, 
which means that for this type of change one can- 
not maintain the micro scale ratios. 

These modifications of physical properties are 
equivalent to changing the constants of nature q 
and c; the unit of charge and the speed of light. 
Then determines the ratio between micro and 
macro scales through the charge density and c^^ 
influences the ratio of the speeds. fiQ needs to 



be kept constant, in order to keep the magnetic 
pressure and hence the plasma beta unchanged, 
hence eo is the parameter that needs to be mod- 
ified. Here is a summary of these two types of 
modifications: 

A: PHYSICAL CHANGES: E - 6, p - 6^, T = 
constant 

CONSEQUENCES: A^ - dskin ^ Vg ^ 
b^^,ujp ^ u]g ^ b,VA ^ Vs ^ vth — constant 

B: PHYSICAL CHANGES: T - t^,g - 
, B ^ t, p = constant 
CONSEQUENCES: A^ - t, but Ssk^n ~ 
ujp — constant, va ^ Vs ^ vth ~ t 

The idea of modifying the cons tants of na- 
ture h as previously been employed. IPrake et al 



()2006l ) r educed in their approa c h the speed of light 
c, while ISiverskv fc Zharkoval (l2009f ) reduced the 
particle density, which has the same effect on the 
micro-scales as lowering the charge per particle. 
In fact, one can show that changes of on the one 
hand temperature and particle density, and on the 
other hand changes of the speed of light and the 
charge per particle, have similar effects: Decreas- 
ing the speed of light or increasing the tempera- 
ture (squared) both reduce the ratio of the speed 
of light to the thermal speed, and decreasing the 
particle density or the charge per particle both 
increase the ratio of micro- to macro-scales (ad- 
justing at the same time the magnetic field so as 
to keep the plasma beta the same). So, in that 
sense, our changes of the speed of light and the 
electric charge per particle effectively corresponds 
to simulating a coronal region with a very high 
temperature and a very low particle density. 

Apart from the modifications discussed above, 
a reduction of the ratio of the electron to pro- 
ton mass from 1836 to 18 is used. This is done 
to decrease the gap between the ion and electron 
plasma frequency and skin depths to acceptable 
values. It is a standard trick in PIC simulations, 
and mass ratios above 16 are normally considered 
enough to separate the two scales. 

Decreasing the ratio between the different 
speeds (e.g. by lowering the speed of light) are 
generally motivated by an assumption that if the 
speeds are nevertheless much smaller than the 
speed of light, one expects only very marginal 
changes in the dynamics, while the savings in com- 
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puting time (which scales as the speed of hght) are 
considerable. 

Changing the charge per particle with a large 
amount (on the order of 10^ here) is more dras- 
tic, but is necessary to bring micro scales into 
a realm that is resolvable with current comput- 
ing resources. To lowest order, such a change of 
scales is not expected to change the magnetically 
dominated dynamics dramatically; charged parti- 
cles are still forced to move essentially along mag- 
netic field lines, with gyration orbits oriented in 
essentially the same manner with respect to the 
large scale structures. What changes is exactly 
what is required to change, namely the ratios be- 
tween micro and macro scales. Any effect that de- 
pends on this ratio then changes in an in principle 
predictable fashion, and one can, a posteriori, at- 
tempt to compensate for this, when analyzing and 
discussing the results, as is done below in Section 



A crucial point when modifying the constants 
of nature is to ensure that the hierarchy of char- 
acteristic speeds, times and length scales are, to 
the largest extent possible, kept as in the real av- 
erage coronal environment. For the characteristic 
speeds this hierarchy is 



Vd < Vth,p ^ Vs < Vth,e <VA < C, 



(8) 



where Vd is the average electron drift speed that 
arises due to the electric current, v^f^ p^^^ is the 
thermal speed of protons(electrons), Vg is the 
speed of sound, va is the Alfven speed, and c 
is the speed of light. For the corona the values 
of these inequalities can be approximated by the 
following numbers [km s~^]: 



0.002 < 90 < 3 900 < 12 000 < 300 000 



(9) 



We performed several simulation runs with differ- 
ent physical resolutions and modifications of q and 
c, while measuring everything in scaled SI units 
(unit of length = 1 km, unit of time = 0.1 s, and 
unit of mass density = 10^^^ kgm^"^). 

In Table [T] the typical initial run parameters for 
a subset of the different simulations performed are 
summarized. The infiuence of the modifications of 
the physical properties in different simulations on 
our results are investigated in Section |4?2] and |4?5] 
In contrast to MHD simulations of a solar event, 
the magnitude of the parameters of a scaled PIC 



simulation cannot be directly compared to obser- 
vationally obtained values, since the modifications 
of the constants of nature have an influence on 
the magnitudes of the parameters. But estimat- 
ing what the corresponding solar values would be 
is possible, by careful rescaling. 

4. Results and Discussion 

There are several objectives with these experi- 
ments, but the most important one is to establish 
whether any particle acceleration (i.e., production 
of non-thermal particles) takes place in the exper- 
iments and if so, what the main particle acceler- 
ation mechanism at work is. As we demonstrate 
and discuss below, non-thermal particles are in- 
deed being produced. Their energy distributions 
are approximate power-laws, with slopes similar to 
the slopes inferred from observations of the bulk of 
sub-relativistic accelerated electrons in magnetic 
reconnection events in the solar corona. 

Most importantly the presented simulations 
manifest the main electron accelerator as being 
a systematic electric field, which evolves relatively 
slowly, and which we therefore characterize as be- 
ing essentially a 'direct current' (DC) electric field. 
This does not mean that the field is completely 
stationary, nor that time evolution is unimportant. 
What it does mean is that we can demonstrate 
that power-law distributions of electrons can be 
created even by nearly stationary electric fields, as 
a result of basically geometric factors, and with- 
out reliance on a recursive process in time. Such 
systematic electric fields occur in our experiments 
in connection with the strong current sheets that 
build up as a consequence of the imposed stress 
on the magnetic field. 

In the beginning of the simulations we observe 
a strong increase in the electric current density 
and the electric field in the fan-plane, close to the 
location where the relative inclinations between 
the magnetic field lines of the two conductivity 
domains — the inner and outer fan-plane area — are 
largest, but also close to the inner spine. Due 
to the boundary conditions chosen we see a mo- 
tion of the inner spine as well as a relative motion 
between the inner and outer spines, where their 
relative distance increases with increasing stress 
and decreases as the system relaxes, as discussed 

The photospheric 



(12012) 
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Table 1: Summary of simulation runs. 



Run 


As [km] 




At 


rg/As 


Ad /As 


c [km s "'"] 


vj [km s "'"] 


# cells 


# part. 


t [s] 


IS 


70 


5 


0.0072 


0.24 


0.4 


3900 


28 


54 M 


2. IB 


21 


2M 


35 


10 


0.0036 


0.48 


0.8 


3900 


28 


427M 


17B 


17 


3M 


35 


5 


0.0036 


0.2 


0.4 


3900 


32 


427M 


17B 


30 


4M 


35 


5 


0.0013 


0.2 


0.2 


7800 


32 


427M 


17B 


17 


5L 


17.5 


5 


0.0018 


0.1 


0.4 


3900 


22 


3.4 B 


136 B 


12 


3Mt 


35 


5.5 


0.0036 


0.77 


1.3 


3900 


63 


427M 


17B 


12 


2Mt 


35 


11 


0.0036 


1.55 


2.6 


3900 


67 


427M 


17B 


12 



Note. — We use the median values of the entire domain for all ratios of parameters/As, where As is the cell size, <5e = c/uipe 
the skin depth, uipe the plasma frequency, Xo = vth/^pe the Debye length, vth,e/'^ce the gyro radius, uice the electron 

gyro frequency, vj = J/ (qn) is the electron current velocity, and At the time step. The simulated time is given in solar seconds. 
The runs with the attribute T have twice the temperature compared to the other simulations runs. 



boundary motion triggers magnetic reconnection 
at the null-point by adding shear between the in- 
ner and the outer magnetic field lines of the fan 
surface, which causes stress on the system that can 
only be reduced by dissipative processes, includ- 
ing magnetic reconnection. The resulting electric 
field generated during these processes is shown in 
Figure [2j which also shows the location of the non- 
thermal high energetic particles (arrows), located 
in the current sheet of the fan-plane outlined by 
magnetic field lines, coinciding with the location 
of an enhanced electric field and E||B. The semi- 
transparent light blue surface depicts the highest 
electric current density, which is also located in 
the fan-plane. The cut-out around the null-point 
is taken after 8 seconds of evolution in simula- 
tion run3M. A fraction of the high-energetic elec- 
trons move through the null-point up along the 
outer spine. Due to the increasing magnetic field 
strength with increasing distance from the null, 
these particles will essentially follow the magnetic 
field lines, and will, unless refiected by a mag- 
netic mirroring effect, likely impact the solar sur- 
face where the strongly bent outer spine intersects 
with the photosphere. The other (larger) fraction 
of non-thermal particles passes the null-point and 
impacts the photosphere on the north-east side of 
the fan-plane (to the left in Figure . A com- 
parison of this impact region with the observa- 
tions of this event reveals several similarities, as 
discussed in Section 14.41 The electrons in the fan- 
plane experiences a strong acceleration parallel to 
the magnetic field inside the current sheet while 



approaching the null-point. 

4.1. The Occurrence of Electric Fields in 
MHD and PIC Experiments 

A systematic electric current develops in the 
fan-plane as a result of the photospheric driver, 
which introduces a near discontinuity in direc- 
tion of the magnetic field in the two regions with 
different connectivity; outside and inside the fan 
plane. This occurs already in the MHD experi- 
ment, where the electric current J is identically 
equal to V x B (disregarding constant factors), 
and where the corresponding electric field in the 
local reference frame, here called the diffusive elec- 
tric field, is given by an Ohm's law of the type 
E,, = 77 J, where 77 is a numerical resistivity. 

Apart from rapid fluctuations (on plasma fre- 
quency time scales) a similar equality between 
V X B and the electric current J must also hold 
in the PIC experiment. Because the PIC simu- 
lations are collisionless one would perhaps expect 
that the diffusive part of the electric field, initially 
not inherited from the MHD simulation (only the 
advective, — u x B, part is kept in the PIC initial 
condition), would remain negligibly small. The re- 
sults of our experiment shows, however, that this is 
not the case. Instead of remaining small, the dif- 
fusive electric field actually grows in magnitude, 
although it generally remains much lower (after 
taken into account the modifications as per Sec- 
tion |42]) than in the MHD simulation from which 
the PIC simulations where started. 
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Fig. 2. — A random selection of non-thermal particles (arrows with color indicating their energy) in the 
current sheet (semi-transparent light blue), which is located in the fan-plane, together with the magnetic 
field (dark blue) for runSM. The spine of the topological structure stretches into the left upper corner (outer 
spine) and reaches down into the center of the dome-shaped fan- plane (inner spine) . 



On closer inspection the reason for the growth 
of the electric field becomes obvious: The charged 
particles that carry the electric current — this is 
primarily the electrons — cannot go where-ever 
they want by following the slightest whim of a tiny 
electric field (in the frame of reference moving with 
the current). Instead, they are in general forced 
to follow magnetic field lines. But even along field 
lines the charged particle motions are not unhin- 
dered; already modest (local or global) increases 
of the magnetic field can force a charged parti- 
cle to become reflected. In addition, systematic 
or fluctuating cross-field components of the elec- 
tric field may cause charged particles to drift away 
from its initial field line, especially in low magnetic 
field strength regions, and in regions with signifi- 
cant shear of magnetic field lines; i.e., precisely in 
regions with significant net electric currents. 

In particular, it is clear from the very definition 
of the concept that magnetic reconnection causes 
a continuously on-going change of magnetic con- 
nectivity, and that this change occurs exactly in 
the place where the largest electric current needs 
to be maintained. 

The component of the electric field along B, 
E||B, corresponds closely to the diffusive electric 
field in MHD when considering the corona. E||B is 
plotted together with the magnitude of E in Vm~^ 
for runs IS, 2M, 3M, 5L, 2Mt and 3Mr in Figure 
[31 In the snapshot from the MHD run used for 



the PIC simulations, we find diffusive electric field 
values of about 29 V m~^ and a total electric field 
of 1666 V m~^. A comparison of the locations of 
high diffusive electric fields in the MHD with re- 
gions of high E||B in the PIC simulations show, 
that the peak values can be found in both cases 
in the current sheet, while when comparing the 
advective electric fields, it reveals that, unlike in 
MHD case, the PIC advective electric fields peak 
again in the current sheet. In the MHD case this 
quantity is high in most of the domain around the 
null-point, while in the PIC case it is only high 
inside the current sheet. Both differences (mag- 
nitude and peak location of the advective electric 
field) are presumably a consequence of lowering 
the elementary charge, imposed by the modifica- 
tions. The lower the charge, the higher the speeds 
of the (lightest) particles have to be in order to 
maintain the same electric current, set by the mag- 
netic field geometry. Therefore the diffusive as 
well as the advective electric field adapt to higher 
values in the PIC simulations, particularly in the 
current sheet. As previously mentioned, our modi- 
fications of the constants of nature prevent us from 
directly comparing these electric field values to the 
values obtained in the MHD simulations. But, as 
shown in the next section, we are able to make a 
qualified guess as to what the PIC electric fields 
would be in the real solar case, which we further 
compare to the MHD electric fields. While most 
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medium resolution runs (cf. Table [T]) show com- 
parable electric fields (Figure |3]), with differences 
of a factor of 2-3, run IS is under-resolved, and 
clearly disagrees with all other simulation runs. 

Charged particles moving in a realistic (non- 
smooth) electro-magnetic field effectively experi- 
ence a 'resistance', and a non-negligible, system- 
atic electric field is needed to maintain the elec- 
tric current J consistent with V x B, as required 
by the Maxwell equations. Below we demonstrate 
that, at least in our numerical experiment, it is this 
systematic electric field that is mainly responsible 
for the particle acceleration. We then argue that, 
since the same phenomena must occur in the real 
solar case, a similar particle acceleration mecha- 
nism must be at work there. 

4.2. The DC Electric Field as the Particle 
Accelerator 

The systematic electric field which develops 
across the fan-plane peaks where the electric cur- 
rent density is largest. This is where particles 
are mainly being accelerated and where the major 
contribution to the power-law energy distribution 
population presented in Section 14.51 comes from. 
Its origin is the tendency for an imbalance between 
the current density and the curl of the magnetic 
field, which occurs as a consequence of the dissi- 




10^ 1 

2 4 6 8 10 12 14 
Time [sec] 

Fig. 3. — Average magnitude of the largest total 
electric field (solid) and E||B (dashed) calculated 
on a cut-out around the null-point with dimen- 
sions 13 X 21 X 16 Mm over the number of grid 
points covering 1 Mm for runs 2M, 3M and 5L and 
the higher temperature runs 2Mt and SM^. 



pative reconnection processes. As magnetic field 
lines are reconnected, the charged particles that 
flow along them become 'misdirected' and need to 
be replaced by the acceleration of new particles 
that are now, instead of the previous ones, situ- 
ated correctly with respect to the magnetic field 
and its curl. 

Tracing particles with kinetic energies in the 
high-energy tail and interpolating the fields to the 
particle position in time steps of 0.003 seconds 
for 5 solar seconds in run 3M, the direct corre- 
lation between a continuous energy gain for elec- 
trons and a negative electric field component in 
the direction of motion is apparent. These non- 
thermal particles are primarily found in regions 
close to, or inside, the current sheet of the fan- 
plane, where the electric current density is high- 
est. In Figure 21 seven representative non-thermal 
particles that gain energy during approximately 
4 seconds are shown. Similar plots, but for par- 
ticles that loose energy, are shown in Figure [5] 
Electrons starting from the right hand side of the 
current sheet feel first a very diffuse and rapidly 
changing electric field, due to a fragmentation of 
the current sheet, which first develops there. This 
behavior is evident when looking at the green en- 
ergy bump around time 5 seconds in the first panel 
of Figure 01 In such fluctuating regions there is a 
constant competition between the dominance of 
the perpendicular electron movement in-between 
the strong E||B patches and the parallel motion 
inside an E|1B region. But once a particle is inside 
a current filament the velocity is mostly directed 
oppositely to the electric field, and electrons then 
experience a rapid acceleration, due to which the 
perpendicular energy relative to the magnetic field 
becomes negligible and the electrons move almost 
parallel to the magnetic field (cf. the cos(pitch an- 
gles) in Figure [Hand Figure [5]). The gyro radius is 
permanently very small (on the order of 3 - 8 km) , 
meaning that electrons tightly follow the magnetic 
field lines. When particles approach the null-point 
(marked by the cross- hair in the graphics), their 
gyro radius naturally grows, as rg o: v± ■ . At 
the same time the magnetic field is weak enough 
that the electrons are no longer strongly confined 
to the magnetic field lines, and they are instead 
directly accelerated by the electric field, rather in- 
dependent of their orientation relative to the mag- 
netic field. Particles entering the region close to 
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the strongest current sheet may also experience 
a modest growth in gyro radius, as the magnetic 
field there is rather weak. 

Small scale regular oscillations which can es- 
pecially be seen in the E||v and in the cos (pitch 
angle) in the panels of Figures |4] and [5] occur at a 
frequency of about 7.65 s~^, which is of the order 
of the plasma frequency. 

We provide in Figure [3] as well as in Figure S] 
and [5] electric field values which are affected by 
our modifications of the constants of nature. We 
therefore need to address how to transform these 
values back, in order to make them comparable to 
the MHD simulation results as well as to observa- 
tions. 

For this purpose we start out with the force that 
an electron with charge q and mass experiences 
when moving in an electric field i?, assuming that 
the DC electric field is the main particle accelera- 
tor F — q E = nie a. We further neglect the dis- 
placement current since the accelerating electric 
field - ignoring small scale turbulence - is persis- 
tent. At time t and after having covered a distance 
L = \at^ the electron has a velocity v of 

v^at^{2Laf^^('^-^\'^ (10) 
\ me J 

The electric current has to balance the magnetic 
field through Ampere's law 

J = A^„-iVxB«M-i— , (11) 

where AB is the typical change in the magnetic 
field across the current sheet and AL is the thick- 
ness of the current sheet. The current is generated 
by moving the electrons 

J = n.q^n('-^\'\ (12) 
\ rrie J 

where n is the particle density. Using the above 
two expressions for J we obtain an equation for 
the total energy gain of a single particle due to 
the electric field 

Given that this acceleration only happens inside 
the current sheet the expression can be interpreted 
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Fig. 4. — Random particles from the power-law 
energy tail of run 3M, which gain energy, is 
the gyro radius, cos{B,v) the cosine of the pitch 
angle and E||v the electric field component point- 
ing in the direction of the particle velocity. The 
background images in the two panels to the right 
are the sum over electric current density slices. 
The upper panel shows the sum of xy-slices from 
a height of 2.5 Mm to 2.9 Mm illustrated in a) 
(see top of figure), while the lower panel has xz- 
slices averaged over 12.6 Mm to 13.4 Mm in the y- 
direction as shown in b) . Both are taken at t = 6 s 
and raised to the power 0.5 to enhance the visibil- 
ity of the fine structures. The null-point is located 
at [26.5, 13.9, 2.8] Mm, plotted as a cross-hair. 
Note that the particle trajectories are projected 
onto these planes, and that the dynamics changes 
slightly during the time interval. Therefore the 
trajectories are seen at a slight displacement. 

as the maximal acceleration a single electron can 
obtain, if it moves continuously inside the current 
sheet over the distance L. 

To relate Equation [13] to our simulations and 
real observations, we need the scaling of AB and 
the current sheet thickness AL. AL is determined 
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Fig. 5. — Analog to Figured but for random par- 
ticles from the power-law energy tail of run 3M, 
which lose energy. 

by stability constraints, as instabilities diffuse the 
current away, in case it shrinks below a certain 
thickness AL. There have been several studies 
in 2D and 2.5D on the current sheet thickness 
concluding that AL is comparable to the size of 
the diffusi on region , as no ted in laboratory exper- 
iments bv iji et al] (2008 ) , as w ell as in PIC sim- 



ulations bv iHesse et al. ( 200l[ ). But these mag- 



netic geometries are not directly comparable to 
our case, in which most of the diffusion takes place 
in the current sheet rather than around an ideal- 
ized magnetic X-point geometry. In fact, very lit- 
tle is known about the thickness of 3D fan-plane 
current sheets. We therefore cover here the most 
probable cases. We assume the smallest length 
scale over which a coherent large scale plane of 
current can be maintained is either the electron 
gyro radius rg in the magnetic field or the elec- 
tron skin depth. 



AL 
AL 



meVth,e 
qB 

/ me 



1/2 



(gyro radius) (14) 
(skin depth) (15) 



(if AL is instead of the order of the ion gyro radius 
it would be larger than the estimate in Eq.[T4]by a 
factor (TOp/me)^/^, which is about a factor of 4 in 
our case) . If we rewrite AB as a fractional change 
in the magnetic field AB = ebB, and express B 



and AB in terms of the Alfven speed va in the 
plasma 

2 ^ ^ 

Va = = ^ 16 

two elegant expressions emerge for Equation I13[ 
defining the maximal electron energy generated by 
the DC acceleration 

qEL = e%—^, (AL - gyro radius) (17) 



th.e 



qEL^e%EA, 



(AL - skin depth) (18) 



where Ea — ■^mpVA is the "kinetic Alfven en- 
ergy" , which is needed to move information in the 
system, and E, 



th,e 



is the thermal en- 



ergy of the electrons. If the current sheet thickness 
is related to the gyro radius, then the higher the 
temperature the smaller the acceleration of an in- 
dividual particle. A higher temperature is also re- 
flected in a larger gyro radius, but the total energy 
available for acceleration induced by Ampere's law 
is the same, and hence there must be more, but 
lower energy, particles in a thicker current sheet. 
On the other hand, if the current sheet thickness 
is related to the electron skin depth the maximum 
energy should be independent of the temperature. 
Note that the right hand side of Equations [17] and 
[THonly contains macroscopic fluid parameters. 

Equations [TTl and fTHl provide an estimate of the 
electric field and at the same time predict what 
the maximum resulting particle energies would 
be if we were able to run without reducing the 
charge per particle; the charge times the electric 
field magnitude is to lowest order a constant and 
hence our modifications in charge are reflected in 
a qmod times too large electric field, where typi- 
cally qmod ~ 2 X 10^. Equations similar to [17] or 
[T8l are what we expect to be able to test in the 
future. As for now the numerical resolution is not 
sufficient for an experiment to resolve typical gyro 
radii with enough grid cells, while simultaneously 
resolving the large scale plasma. 

In the current experiments the current sheet 
thickness is essentially determined by the grid 
spacing, and is generally a few times As; about an 
order of magnitude larger than the typical electron 
gyro radii, and about half an order of magnitude 
larger than the typical proton gyro radii. If, at 
the same elementary charge per particle q, we were 
able to increase the numerical resolution in order 
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to resolve the gyro radii and the current sheet be- 
came correspondingly thinner, then according to 
Equation ll3l the maximum energy gain would in- 
crease with the square of the AL factor, so with 1 

- 2 orders of magnitude. 

We conclude that a conservative estimate of the 
current sheet electric field in the Sun would be 
smaller by at most the factor (about 2 x 10^) by 
which the elementary charge per particle has been 
reduced, and it could possibly be 1 - 2 orders of 
magnitude larger. Taking 4000 Vm~^ as a typical 
magnitude of the electric field in our experiments, 
we thus come up with a conservative estimate on 
the order of 2 mV m""'^ for the solar electric field in 
a situation analogous to the one we model (which 
is not a flaring situation). Using a less conserva- 
tive estimate based on the argument of numerical 
resolution, the electric field is on the order of 20 

- 200niVm^^. In comparison electric fields dur- 
ing solar flares are i nferred to be on the order of 
thousands ofVm"! (jOiu et al.ll2002l ). 



4.3. 



Magnetic Field Geometry versus Power- 
lav^r Distribution 



As the current sheet channels make up a very 
small fraction of the domain, most particles will 
not exhibit the correct angle to exactly pass 
through an acceleration channel, but will instead 
be deflected by the magnetic field, ending up out- 
side the current sheet without undergoing a con- 
tinued acceleration. Hence most particles are not 
continuously accelerated. The electric current it- 
self is mainly carried by the lowermost part of the 
power-law distribution, as illustrated in Figure [Bl 
which shows the contributions from the particles 
to the electric current density from three equally 
large regions of log(energy) of the power-law tail 
of run 5L. The particles of the lowest bin are most 
probably in constant exchange with the thermal 
particle distribution. Figure [7] sets the power-law 
tail in relation to the bulk flow and the thermal 
particles. The power-law tail dominates over the 
thermal contribution with respect to contribu- 
tions to the electric current density. About 0.2% 
of all electrons in the computational box make up 
the non-thermal electrons high-energy tail of the 
distribution. 5% of the total electron energy is 
carried by the power-law tail particles, while most 
of this energy is in the most energetic particles of 
the power- law tail, as illustrated in Figure [5] But, 



if we only consider particles with a negative ver- 
tical velocity, thus moving towards the bottom of 
the box, and additionally reside in a zone on the 
lower quarter of the box of 1.9x10^^ km^ (from 
approximately 0.175 - 1.750 Mm above the bot- 
tom boundary) , the energy share coming from the 
power-law tail amounts to over 50% and we find 
more than 5% of all particles in the high-energy 
tail population in run 5L, while the total num- 
ber of electrons in this cut-out is about 1.52x10^. 
These are the particles that on the real sun would 
be decelerated in the chromosphere, leaving an 
imprint in the form of observable bremsstrahlung 
emission. 

4.4. Comparison with Observations 

Considering the observations of this particular 
reconnection event, in addition to the bright rib- 
bons observed in Ha and UV at the intersection 



of th e fan and the chromosphere (jMasson et al 



20091) . the chromospheric foot points of the inter- 
change reconnection region show also soft and 
hard X-ray signatures in the area northwards 
of the null-point, presumably originating from 
bremsstrahlung (H. Reid, private communica- 
tion). This coincides well with the impact region 
of the power-law electrons in our simulation which 
travel along the fan-plane and finally hit the lower 
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X [Mm] X [Mm] X [Mm] 

Fig. 6. — Run 5L: Power-law tail electric current 
density (J) contributions from particles within the 
energy range of a) 4.3x10"^^ < e < 1.4x10^^^, b) 
1.4xl0~2i < e < 4.4x10^21 ^nd c) e > 4.4x10^21 
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24 26 28 30 26 28 30 26 28 3C 
X [Mm] X [Mm] X [Mm] 

Fig. 7. — The electric current density arising from 
a) bulk flow b) the power-law tail particles and c) 
the thermal particles in run 5L. 



boundary. Figure [9] shows the non-thermal elec- 
tron energies at their impact regions on the lower 
boundary of the box accumulated over t = 4 - 
9 seconds. The small difference in the location 
of the peak intensit y compared to observations 



( Masson et al. 2009() can at least partly be ex- 



plained by the driving pattern, which is of course 
only an approximation to the real photospheric 
boundary motion. A second reason may be the 
overall magnetic field at the start of the simula- 
tion, which is an outcome of an MHD simulation, 
which was again initialized from a potential field 
extrapolation. 

However, our electron energies are clearly lower 
than what is needed for the observed emission 
spectra to be produced. Hence it is important to 
emphasize that we do not model the observed flare 
event, but rather the pre-flare phase. One rea- 
son is that we use an MHD state taken at a time 
well before the flare event. Another reason is, as 



shown in lBaumann et al.l ()2012l) . that the bound- 
ary driver in the MHD simulations, representing 
the observed horizontal magnetic fleld motion in 
the active region, does not provide enough shear 
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Fig. 8. — Run 5L: Power-law tail energy (e) con- 
tributions from particles within the energy range 
of a) 4.3x10-22 < e < 1.4xl0~2i^ 1.4x10^21 
< e < 4.4x10^21 and c) e > 4.4x10^21 [a.u.]. 




Fig. 9. — Impact area of non-thermal electrons 
from run 5L. Shown is the electron energy (from 
orange increasing to yellow) added up at the lower 
boundary over a period of 5 sec. Additionally mag- 
netic field lines (purple) passing close to the null- 
point are plotted. 
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and stress to the system to result in an abrupt 
energy release. Despite that, we expect the quali- 
tative nature of the acceleration mechanism to be 
the same in the flare phase as in its pre-phase. 
This is supported by the fact that the highly ac- 
celerated flare electrons are observed to impact at 
essentially the same locations as the non-thermal 
particles in the simulation. 

4.5. The Energy Distribution and the In- 
fluence of the Modifications 

Looking at the energy histogram of the same 
particles considered in Section 14.31 downward 
moving in a cut-out of the lower part of the sim- 
ulation box (see top illustration in Figure [T0|) . we 
find a Maxwell-Boltzmann distribution combined 
with a dN/d\nE = EdN/dE power-law index of 
about -0.78, corresponding to a dN/dE distribu- 
tion power-law index of -1.78 (see Figure ITU)) . A 
power-law index of -1.78 implies that the electric 
current resulting from the power-law population is 
mainly carried by the low energy electrons, while 
the kinetic energy is mainly carried by the high- 
est energy constituents, visualized in Figure [5] 
The color code in Figure [10] shows the temporal 
evolution of the energy distribution function for 
downward-moving particles in the cut-out. The 
power-law index for the full simulation box is sim- 
ilar. Figure [TU] further shows that the tail slope 
rapidly converges toward the power-law index of 
about -1.78, indicating an impulsive acceleration 
of electrons. The weak, apparently non-thermal 
tail present in the distribution function from the 
very beginning (dark blue line), does not share 
origin with the power-law tail that is created dy- 
namically at later times. It arises due to the 
electric current in the MHD current sheet (the ini- 
tial electron velocities are drawn randomly from a 
Maxwellian distribution, shifted with a systematic 
velocity to maintain the local electric current den- 
sity). It is entirely due to the high values of the 
drift velocity necessitated by our rescaled units, 
but does not influence the later creation of accel- 
erated particles; they appear also in tests where 
the electrons are not given a systematic initial 
velocity. 

A comparison of the different electron energy 
distributions from the simulation runs, listed in 
Table [1] may be found in Figure [Tl] and Figure [121 
For these plots all electrons in the full simulation 
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Fig. 10. — Energy distribution of run 5L covering 
12 solar sees (from dark blue to red) and being 
obtained from 1.52x10^ electrons in a volume of 
1.9xl0^^km'^ shown in the illustration above the 
plot. At is 1 solar sec. The dashed line represents 
a Maxwell-Boltzmann distribution fit. 
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Fig. 11. — Energy histogram for all electrons in 
the full box for different simulation runs at time 
t = 12 sec. 
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Fig. 12. — Energy histogram for all electrons in 
the full box at time t = 12 sec for the simulation 
runs with double the temperature of the runs pre- 
sented in Figure [TTl 



domain are used. The histograms are not normal- 
ized by the number of particles in order to allow 
for a better power-law index comparison between 
the different runs. The thermal distribution is the 
same for all simulations having initially the same 
temperature profile and it is stable over at least 
12 solar seconds, as shown as an example in run 
5L, Figure [TOl The non-thermal energy tail part 
quickly approaches a power-law index of about - 
1.75 for all runs, independently of their resolution. 
The power- law itself is mainly a consequence of the 
DC systematic electric field particle acceleration, 
while the cutoff of the energy histogram presum- 
ably is a result of the maximum available electric 
potential difference in the system, determined by 
the current sheet thickness, which is controlled by 
the magnetic field geometry and its evolution. 

For simulations with initially a higher temper- 
ature (see Figure [T2|) the non-thermal part of the 
energy distribution drowns in the thermal part. 
With a plasma beta in the corona << 1, a change 
in the temperature is not expected to significantly 
change the overall electric field (also see Figure 
131) and thereby the maximum energy that can be 
gained by a particle. 

Comparing our power-law index to observations 
shows that our electron power-law index is slightly 
harder than what is indicated by observations, 
where it is often > 2. There have been a num- 
ber of t est particle inves t igations of current sheets, 
such as Turkmani et al. ( 20061) : Wood fc Neukirch 



(|2005l ): IZharkova fc Gordovskvvl (|2005l ) finding 
similar electron power-law indices - and even 
slightly harder. But such comparisons are dan- 
gerous, as the these studies were conducted using 
the test particle approach and the latter addition- 
ally assumed a simplified 3D magnetic and electric 
field configuration, presumably having a signifi- 
cant infiue nce on the power-law index as p artially 
studied bv lZharkova fc Gordovsk^ (|2005t) . 



5. Conclusions 

In this study of a prc-fiarc 3D reconncction re- 
gion near a 3D magnetic null point in the solar 
corona the main electron acceleration mechanism 
has been shown to be the parallel electric field, 
building up as a consequence of the magnetic re- 
connection and dissipation close to the null-point 
and in the fan-plane current sheet. We estimate 
from the simulation results the average electric 
field in this quiescent region to be on the or- 
der of 20 - 200 mV m~^ in the Sun, which is 
significantly smaller (2-3 orders of magnitude) 
than what has been f ound in MHD simulations 
( Baumann et al. 2012f ). but still sufficient to cre- 
ate a noticeable non-thermal electron power-law 
tail. We also present a discussion of the expected 
dependency of the electric field on the modified 
elementary charge per particle; showing that the 
product E ■ q is expected to be conserved to lowest 
order. The verification of this relation must be left 
for future studies due to the current restrictions on 
computational resources. 

The different modification options given in this 
article serve to reduce the ratio between micro and 
macro scales. The main parameter change used to 
increase micro scales is a reduction of the elemen- 
tary charge per particle, essentially equivalent to 
a reduction of the particle density. 

From studies of the locations of the non-thermal 
electrons and of their acceleration paths we con- 
clude that the magnetic field geometry and its 
temporal evolution are likely to be the main fac- 
tors controlling the power-law index measured in 
this experiment. We further show that in the 
lower part of the computational box the electron 
energy is predominantly in the non-thermal com- 
ponent, with a particle impact area that correlates 
well with the observations by TRACE, SOHO and 
RHESSI. In spite of the modifications of the con- 
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stants of nature necessary to make this experiment 
possible we expect that the physical processes in 
the experiment are qualitatively similar to those 
in the real Sun, and we therefore see these kinds 
of studies as a very valuable tool for studying the 
coupling between kinetic and MHD scales in semi- 
realistic models. 

In the future wc hope to strengthen our find- 
ings with even higher resolution simulations, and 
we additionally plan to study the temporal depen- 
dence of the power-law index by employing longer 
simulation runs, which would at the same time 
provide an opportunity to study ion acceleration. 
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